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Roles of orbital degree of freedom of Mn ions in phase transition as a function of temperature and 
hole concentration in perovskite manganites are studied. It is shown that the orbital order-disorder 
transition is of the first order in the wide region of hole concentration and the Neel temperature 
for the anisotropic spin ordering, such as the layer-type antiferromagnetic one, is lower than the 
orbital ordering temperature due to the anisotropy in the orbital space. The calculated results 
of the temperature dependence of the spin and orbital order parameters explain a variety of the 
experiments observed in manganites. 
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I. INTRODUCTION 

Perovskite manganites J 4i_ I 5 J) Mn03 (A: La, Nd, Pr, 
B: Sr, Ca) and their related compounds have been re- 
cently studied extensively from both experimental and 
theoretical sides. Anomalous electric, magnetic and 
structural properties accompanied with the phase tran- 
sition such as colossal magnetoresistance (CMR) attract 
much attention.ETQ Gigantic decrease of the electrical re- 
sistivity is brought about at the vicinity of the transi- 
tion from spin, charge and orbital ordered phase to ferro- 
magnetic metallic one with slightly changing temperature 
and/or applying external fields. ErB It is now accepted that 
such dramatic phenomena are not understood within the 
simple double exchange scenario.!] 

The orbital degree of freedom in Mn ions is one of 
the convincing candidates to bring about not only a rich 
variety of phenomena but also dramatic ones. Due to 
the strong Hund coupling and the crystalline field, two 
e g orbitals of a Mn ion, i.e. the d^ z 2_ r 2 and d x 2_ y 2 
orbitals, are degenerate and one of them is occupied 
by an electron in a Mn 3+ ion. It is well known that 
the (d 3a .2_ r 2 /d 3v a_,.2)-type orbital ordered state, where 
the two orbitals are alternately aligned, is realized in 
the undoped manganites LaMn03.Bil3 After the dis- 
covery of CMR, much study of the orbital states in 
doped manganites has been done. For example, the 
layer (A)-type antiferromagnetic (AF) metal accompa- 
nied with the uniform alignment of the d x 2_ y 2 .orbital is 
observed in Pro. 5 Sro.5Mn0 3 , NaWSr . 55 MnO 3 EHi£l and 
Lai-^Sr^MnOs with x ~ 0.55.t3 Nevertheless, roles of 
the orbital in phase transition are still far from our un- 
derstanding. Since CMR appears near the orbital order- 
disorder transition, it is indispensable to study nature 
of the phase transition in doped manganites where the 
orbital degree of freedom is taken into account. 

In this paper, we study the phase transition in per- 
ovskite manganites based on the model where the or- 
bital degree of freedom and electron correlation are in- 
cluded. By adopting the mean field theory, the phase 
transitions are investigated as a functions of tempera- 
ture (T) and carrier concentration (x). Since there is a 
strong anisotropy in the orbital space unlike the spin one, 



the orbital order-disorder transition is of the first order 
in the wide range of x, and the Neel temperature TV for 
the anisotropic spin ordering, such as the A-type AF one, 
is lower than the orbital ordering temperature Too ■ The 
calculated temperature dependence of the spin and or- 
bital order parameters explains a variety of experiments 
in several manganites. 

In Sect. II, the model Hamiltonian is derived and the 
mean field theory at finite T is introduced. In Sect. Ill, 
the numerical results of the spin and orbital phase dia- 
gram are presented. We focus on the phase transitions 
in (i) lightly doped region where the ferromagnetic Curie 
temperature, Tq, and Too are close with each other and 
(ii) highly doped one where the A-type AF state accom- 
panied with the d x 2_ y 2 orbital appears. In Sect. IV, 
phase transition is studied analytically by expanding the 
free energy with respect to the spin and orbital order 
parameters. Section V is devoted to the summary and 
discussion. 



II. MODEL AND FORMULATION 

Let us set up the model Hamiltonian describing the 
electronic structure in manganites. We consider the 
tight-binding Hamiltonian in the cubic lattice consisting 
of Mn ions. At each site, two e g orbitals are introduced 
and tig electrons are treated as a localized spin (5* 2 ) 
with S — 3/2. We introduce three kinds of Coulomb in- 
teraction between e g electrons at the same site, i.e. the 
intra-orbital Coulomb interaction (U), the inter-orbital 
one (U 1 ) and the exchange interaction (I). The hop- 
ping integral between site i with orbital 7 and its nearest 
neighboring site j with 7' is denoted by t]^ . The Hund 
coupling (Jh) between e g and ti g spins and the antiferro- 
magnetic superexchange (SE) interaction (Jaf) between 
nearest neighboring ti g spins are also introduced. Among 
these parameters, the intra-site Coulomb interactions are 
the largest. Thus, by excluding the doubly occupied state 
in the e g orbitals, we derive the effective Hamiltonian de- 
scribing the low energy electronic states:E3 

H = Ht + Hj + Hh + Ti-AF- (1) 
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The first and second terms correspond to the so-called fr- 
aud J-terms in the t J -model and are given by 
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respectively. Here, t/ = cos (^-mi)T iz — sin( ^m;)T ix 
and (m x ,m y ,m z ) = (1,-1,0). / denotes a direction 
of the bond connecting site i and site j. d ila is the 
annihilation operator of the e g electron at site i with 
spin a and orbital 7 with excluding double occupancy of 
electron and Hi is the number operator defined as rij = 

J~)~„ dL-di-vcr- The explicit form of £77 is determined 
by the Slater-Koster formulasJla 5, is the spin operator 
of the e g electron with 5 = 1/2 and Tj is the pseudo- 
spin one for the orbital degree of freedom defined as 
Ti = (1/2) E 77 ' CT ^7-^)77' 4 where T l2 = +R1/2 
corresponds to the state where the d 3z 2_ r 2 (d x 2_ y 2) or- 
bital is occupied by an electron. Ji = t\j{XJ' — I) and 
,/ 2 = t%/(U' + I + 2Jh) where to is the hopping inte- 
gral between nearest neighboring d^ z 2_ r 2 orbitals in the 
z direction and U = U' + I is assumed. The third and 
fourth terms in Eq. ([!]) represent the Hund coupling and 
the antiferromagnetic SE interaction, respectively, and 
are given by 
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The detailed derivation of the Hamiltonian is presented 
in Ref. [l7[ Characteristics of this Hamiltonian are sum- 
marized as follows; 1) the two kinds of magnetic inter- 
actions between spins of e g electrons, i.e. the SE and 
double exchange-interactions are described by TLj and 
Tit, respectively^ 2) there is a strong anisotropy in the 
pseudo-spin space unlike the spin space, and 3) the first 
term in Eq. (|^) is the dominant term in Tij and stabilizes 
the ferromagnetic state associated with the antiferro-type 
orbital ordered one-where different types of orbital are al- 
ternately aligned .ETa 

Being based on the Hamiltonian, the spin and orbital 
states are studied at finite,-^ and x. The mean field the- 
ory proposed by de Gennescil is applied to the present sy.Sc. 
tern with orbital degeneracy and electron correlation.Lj 
In this theory, the spin and pseudo-spin operators are 
treated as classical vectors: 
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where 9^ is the polar angle in the spin (orbital) space 
and 4>l is the azimuthal one in the spin space. A motion 
of the pseudo-spin is assumed to be confined in the xz 
plane and 6\ describes the orbital state at site i as follows 

I^HcosllS^-O+sinl^-y 2 ). (7) 

From now on, Si and Ti are denoted by vii in the uniform 
fashion. A thermal distribution of Ui is described by the 
distribution function: 



<( K i ) = i r exp(A?.^ 
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where A" is the mean field and vf is the normalization 
factor given by 
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The free energy is given by a summation of the expecta- 
tion value of the Hamiltonian and entropy: 



T= (H) st -TN{S S +S l 



(11) 



where N is the number of Mn ions and S u is the entropy 
for u defined by 



-(\nw u (u)) l 
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(A) u implies the expectation value of A with respect to 
the distribution function. M u = (u ■ X u / (\u\\X u \)) u is 
adopted as an order parameter in (Tij) s t, (Hn)st an d 
(HAF)st- The relation 3(Si) s = (St 2g i) s is assumed. As 
for (Tit) st, the rotating frames in the spin and pseudo- 
spin spaces are introduced. The electron annihilation 
operator is decomposed as di 1<y = z^z* /ij where h\ is 
the creation operator of a spin- and orbital-less fermion 
describing a hole motion and z slyt ) * is an element of 
the unitary matrix in the spin (pseudo-spin) frameEj'EZI 
These are given by z^ — cos(6*f /2)e -i< ^/ 2 , zf^ — 

sin(6' l s /2)e^'/ 2 , z\ a = cos(6»*/2) and zj b = sin(0f/2). H t 
is rewritten as 

H t = (zt^izpV'z^hiht+H.c,, (13) 
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FIG. 1. The spin and orbital phase diagram as a function of 
hole concentration (x) and temperature (T). The schematic 
orbital states are also shown. Tc, Tn and Too indicate fer- 
romagnetic Curie temparature, Neel temperature and orbital 
ordering temperature, respectively. F, A and G indicate fer- 
romagnetic, A-type antiferromagnetic and G-type antiferro- 
magnetic phases, respeceively. The solid (broken) lines is 
for the second (first) order phase transition. The parame- 
ter values are chosen to be Ji/to ~ 0.25, Ji/to = 0.075, and 
JAF/h = 0.0035. 

and is diagonalized in the momentum space as follows 
/ N > \ 

= (EE4/H4-^)> > (i4) 
x k i=i 1 st 

where eL is the energy in the i-th band for the spin- 

and orbital-less fermion with momentum k and Ni is the 
number of the band. The Fermi energy ef is determined 
by the condition x — (1/N) J2k /^(^S ~ £ f) where 
/i?(e) is the Fermi distribution function. The mean field 
solutions are obtained by minimizing T with respect to 

III. NUMERICAL RESULTS 

A. spin and orbital phase diagram at finite 
temperature 

The spin and orbital states at finite T and x are nu- 
merically calculated by utilizing the mean field theory in- 
troduced in the previous section. The four types of spin 
structure, that is, the ferromagnetic (F) structure and 
layer (A)-type, rod (C)-type and NaCl (G)-type AF ones 
are considered. As for the orbital, the ferromagnetic-like 
structure where one kind of orbital exists and the G-type 
one where two kinds of orbital are alternately aligned 
in the [111] direction are considered. The numerical re- 
sults are shown in Fig. ^ Parameter values are chosen 
to be Jx/to = 0.25, J 2 /t = 0.075 and JafIU = 0.0035. 



The schematic pictures of the orbital ordered states are 
also presented. With increasing x from x = 0, the spin 
structure changes as A-type AF — > F — > A-type AF — > 
G-type AF which is associated with change of the or- 
bital states; in the region of x < 0.25, the interaction 
caused by TLj is the dominant one between nearest neig- 
boring orbitals and the G-type AF orbital ordered state 
is brought about. The type of the orbital is denoted as 
(eyey = (f/ - |) where Q\ (B) is the angle of the 
pseudo-spin in sublattice A(B) and its definition is the 
same with 6\ in Eq. (|7|). These orbitals are mixtures 
of d 3x 2_ r 2 and d y 2_ z 2, and d 3y 2_ r 2 and d z 2_ x 2, respec- 
tively. Above x = 0.25, the F-type orbital ordered state 
is realized. In particular, the d x 2_ y v orbital is uniformly 
aligned in the A-type AF spin phase above x — 0.6. A 
large hopping integral for electrons in the xy plane in 
this orbital ordered state is energetically favored in the 
A-type AF state where the hopping in the z direction 
is prohibited.E3 The calculated results of spin and or- 
bital phase diagram at T = are consistent with those 
obtained by the Hartree-Fock theory.Ej The spin and or- 
bital phase diagram at finite T was calculated in Ref. |2^ 
where the Monte Carlo method in a finite-size cluster was 
used in the spin-orbital-lattice coupled model. 

Now, let us focus on the spin ordering temperatures, 
i.e. Tc and Tjv, and the orbital ordering one, Too, 
in Fig. |l|. These ordering temperature vs. x curves 
qualitatively reproduce the,-£xperimental res-ultsj ob- 
served in Lai_^Mn0 3 ,EM3 Pri-^Sr^MnOaBlla and 
Ndi-zSr-jMnOstj except for the narrow region of the 
charge ordered phase in Ndo.sSro.sMnOa. It is shown in 
Fig. 1 that Too is higher (lower) than Tc in the region 
of x < 0.1 (x > 0.1). The dominant interaction in the 
region of a; < 0.1 is provided by 7ij where the effective 
interaction between orbitals in the paramagnetic state 
and that between spins in the orbital disordered state 
are given by 3Ji/2 and Ji/2, respectively. Here, the first 
term in Hj is considered. Thus, Too is higher than Tc- 
On the other hand, in the region where Tit is dominant, 
gain of the kinetic energy associated with the long range 
ordering causes the transition; it is assumed that doped 
holes are introduced at the bottom of the band and the 
kinetic energy is proportional to the band width. The ra- 
tio of the band width in the ferromagnetic state to that 
in the paramagnetic one is obtained as 3/2 where the 
orbital disordered state is assumed. On the other hand, 
the ratio of the band width in the F-type orbital ordered 
state to that in the orbital disordered state is obtained 
as 7r 2 /8 where the spin paramagnetic state is assumed. 
The energy gain associated with the orbital ordering is 
smaller than that with the spin one, so that Tc is higher 
than Too- This is attributed to the hopping integral 
between different kinds of orbital. 

Between Tjv for the A-type AF state and Too , the re- 
lation Tn < Too is satisfied in the whole region of x in 
Fig. [l| In addition, the orbital order-disorder transition 
and the A-type AF one are of the first order in the region 
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FIG. 2. (a) Temperature dependence of the magnetization 
Mf (solid line) and the orbital order parameter Mt (broken 
line) at x — 0.14. The inset shows the magnetization curve 
in Lao.88Sro.12MnO3.Ej (b) Magnetic field dependence of the 
orbital ordering temperature Too at x = 0.14. The inset 
shows Too at gfisB/to — and 0.01 and the ferromagnetic 
Curie temperature Tc around x c - Parameter values in the 
calculation are the same as those in Fig. |lj 



of x > 0.25. This is numerically confirmed by disconti- 
nuity in the orbital (spin) order parameter at Too CT/v)- 
Both the two results originate from the anisotropy in the 
pseudo-spin space as discussed latter in more detail. 



neighboring spins is given by — 2 Jx((niTij)/4 — ( T iT l j)) 
where the first term in Ttj is considered. With taking 
into account the fact that the orbital ordered state is 
(6a, ©b) = ( 7r /,37r/2) in this region of x, the effective 
interaction is rewritten as —2J\{1 — x) 2 (l/4 + 3M 2 /16) 



for I = x and y and — 2Ji(l — x) /4 for I = 



The 



effective magnetic interaction and thus the the magne- 
tization increase below Too- I n the opposite way, the 
orbital order-disorder transition is influenced by change 
of the spin state. The magnetic field dependence of Too 
at x — 0.14 is presented in Fig. 2 (b). The inset shows 
change of the phase diagram by applying the magnetic 
field around x c . Too increases with applying a magnetic 
field. The effective interaction between nearest neighbor- 
ing orbitals is given by 2J 1 (1 - x) 2 (3/4 + M 2 /4) which 
increases by applying the magnetic field. This implies 
that the orbital state is controlled by the magnetic field, 
although the field is not a canonical external field for the 
pseudo-spin. 

This unique phase transition originating from the 
coupling b etwe en spin and orbital is observed in 
manganitesJljH In Lao.88Sro.i2Mn03, the ferromagnetic 
ordering occurs at 175K and the orbital ordering is con- 
firmed below 145K by the resonant x-ray scattering which 
is a direct probe to detect the orbital ordering. The ferro- 
magnetic phase with the orbital disordered state changes 
into the phase with the orbital ordered state at 145K, 
so that it corresponds to the calculated phase transition 
at Too around x = 0.14. It is experimentally confirmed 
that the magnetization is enhanced below 145K (the in- 

srature 
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set of Fig. 2(a)) and the orbital orderin; 
increases with applying the magnetic field 
experimental results are well explained by the present cal- 
culation and are strong evidences of the novel coupling 
between spin and orbital in this compound. 



B. spin and orbital phase transitions in lightly 
doped region 

In the lightly doped region in Fig. 1, Tc (Too) in- 
creases (decreases) with increasing x from x = and 
the two transition temperatures cross with each other 
around x = 0.13 termed x c . Around x c , the coupling be- 
tween spin and orbital degrees of freedom brings about 
the unique phase transition as follows. Let us focus on 
the region where x is slightly higher than x c . With de- 
creasing T, the system changes from the paramagnetic 
phase with the orbital disordered state to the ferromag- 
netic phase and then to the ferromagnetic one with the 
orbital ordered state. The temperature dependence of 
the magnetization at x = 0.14 is presented in Fig. 2 (a) 
where the order parameter of the orbital ordered state is 
also plotted. It is shown that the magnetization is en- 
hanced below Too- This originates from the coupling 
between spin and orbital in Jij (Eq. (|^)); in the mean 
fields theory, the effective interaction between nearest 



C. phase transition in highly doped region and 
A-type AF metal 

In this subsection, we focus on the phase transition in 
the highly doped region (x > 0.5) in Fig. 1 where the A- 
type AF state associated with the d x 2_ y 2 orbital ordered 
one appears. There are two kinds of carrier concentration 
regions termed region I (0.8 > x > 0.6) where a sequen- 
tial phase transition from the paramagnetic state to the 
ferromagnetic one and to the A-type AF one occurs with 
decreasing T, and region II (x > 0.8) where the transi- 
tion from the paramagnetic state to the A-type AF one 
occurs. The temperature dependences of the spin order 
parameters at x = 0.725 (region I) and x = 0.9 (region 
II) are shown in Figs. 3(a) and (b), respectively. Mf 
and Mf F are the order parameters of the ferromagnetic 
and A-type AF spin structures, respectively. As shown 
in Fig. 3(a), Mf appears at Tc/t — 0.95 where the 
transition is of the second order. With decreasing T, the 
F-type orbital ordered state with d x 2_ y 2 orbital appears 
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FIG. 3. Temperature dependences of the spin order pa- 
rameters (a) at x — 0.725 and (b) at x = 0.9. The solid and 
broken lines show the order parameters in the ferromagnetic 
state {Ms) and the A- type antiferromagnetic one (M^ F ), re- 
spectively. Parameter values are the same as those in Fig. 
[l The inset in (a) shows the temperature dependences of 
magnetic Bragg reflections in Pro.5Sro.5MnO3.tj 



a t Too Ao = 0.72 and the transition from the ferromag- 
netic phase to the A-type AF one occurs at Tn- This 
transition is of the first order and the canted AF phase 
does not appear between the ferromagnetic and A-type 
AF phases. This sequential phase transition is caused by 
the thermal fluctuation of the orbital; as previously men- 
tioned, the A-type AF state and the F-type orbital or- 
dered one with d x 2_ y2 orbital are cooperatively stabilized 
at T = 0. With increasing T, the thermal fluctuation of 
orbital grows up and the hopping integral in the z di- 
rection becomes finite. As a result, the double exchange 
interaction in this direction overcomes the antiferromag- 
netic SE one and the ferromagnetic phase is stabilized. 
In the region II, the A-type AF ordering and the F-type 
orbital one with d x 2_ y 2 orbital occur at the same temper- 
ature where the transition is of the first order as shown 
in Fig. 3(b). In both the regions I and II, the relation 
Tn < Too is satisfied. The first order transition at Too 
and this relation between Too and Tn originate from the 
breaking of the inversion symmetry of the system with re- 
spect to the orbital pseudo-spin operator, as discussed in 
the next section. The calculated results of the sequential 
phase transition in the region I well reproduce, the exper- 
imental results observed in Pro.sSro.sMnOstS (the inset 
of Fig. H (a)). The first order transition at T o{= T N ) 
in the region II is consistent with the experiments in 
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FIG. 4. Temperature dependences of the spin and orbital 
order parameters (a) at x = 0.725 and (b) at x = 0.9 in the 
tetragonal lattice. The solid, broken and dotted lines show 
the order parameters for the ferromagnetic structure (Mf), 
the A-type AF one (M^ F ) and the F-type orbital ordered 
one {Mt), respectively. The uniaxial anisotropy is introduced 



as £7*5 = y^JTM 



R with R = 1.2 where t 



xy(z) 



and 



Jaf are t ne hopping integral and the antiferromagnetic SE 
interaction in the xy plane (the z direction). The other pa- 
rameter values are the same as those in Fig. |l[ 



Nd .45Sro.55Mn0 3 BEj where the Mn-0 band length in 
the z direction (xy plane) is confirmed to become short 
I It implies the F-type orbital ordering with 



(long) at T N £ 

d x 2_ y 2 orbital as predicted from the present calculation. 

In the actual compounds, where the A-type AF state 
is observed, the tetragonal lattice distortion is observed 
and the cubic symmetry is broken far above Tn oo We 
simulate this distortion by introducing the uniaxial anis- 
toropy of the hopping integral and the SE interaction 
and investigate the phase transition. It is assumed that 

*o7*o = VJ'af/Jaf = R wher e C W and J^ z) are 
the hopping integral and the antiferromagnetic SE inter- 
action in the xy plane (z direction). The temperature 
dependences of the spin and orbital order parameters at 
x = 0.725 and x = 0.9 with R — 1.2 are shown in Figs. |] 
(a) and (b), respectively. Due to the uniaxial anisotropy, 
M t is finite above the orbital ordering temperature in 
the system with the cubic symmetry. It is worth to note 
that 1) Tjv for the A-type AF state increases and 2) the 
transition at Tn is of the first order, although the dis- 
continuity of Mf F is reduced. The latter is attributed 
to diminution of the change of M t at Tn- 

Finally, the magnetic field dependence of Too is shown 
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FIG. 5. Magnetic field dependence of the orbital ordering 
temperature at x — 0.9 and schematic spin configurations. 
Parameter values are the same as those in Fig. 



in Fig. 5. Too decreases with applying the magnetic field 
in the region of g/iB-B/fo < 0.0025 where the nearest 
neighboring spins in the z direction are canted. The spins 
become parallel at g/igB/ta = 0.0025 termed B c , and 
Too increases with increasing the magnetic field above 
B c . The orbital ordered state below Too is of the F-type 
with d x 2_ y 2 orbital and does not depend on magnitude 
of the magnetic field. Above and below B c , the differ- 
ent mechanisms dominate the magnetic field dependence 
of Too', in the region of B > £> c , the spin canting due 
to the magnetic field promotes the electron hopping in 
the z axis which weakens the orbital ordered state. On 
the other hand, above B c , magnitude of the magnetic 
moment is enhanced by increasing of the magnetic field 
and the orbital ordered state associated with the ferro- 
magnetic spin structure is stabilized. The magnetic field 
dependence of T/v for the A-type AF spin structure is re- 
cently measured in Ndo.45Sro.55MnO3.Eil It is experimen- 
tally shown that T/v gradually decreases with increasing 
the magnetic field. From the present calculation, it is 
predicted that this reduction of T/v is accompanied with 
that of Too, and with increasing the magnetic field fur- 
thermore, Too increases above the critical value of the 
field. 



IV. FIRST ORDER TRANSITION AND ORBITAL 
DEGREE OF FREEDOM 

As mentioned in Sect. Ill, several characteristics of the 
phase transitions in manganites are attributed to the 
unique properties of the orbital degree of freedom. In 
this section, we study analytically the phase transition 
by expanding the free energy— with respect to the spin 
and orbital order parameters.! 2 .! Let us consider the fer- 
romagnetic and A-type AF spin structures and the F 
and G-type AF orbital ones. The expectation values of 
the Hamiltonian at finite temperature are calculated by 



the mean field theory introduced in Sect. II and are ex- 
panded with respect to M s and M t . It is assumed that 
doped holes are introduced at the bottom of the band 



denoted by e~ in Eq. (|14|). The explicit form of (Tij} s t is 
given by 

(Hj) st _ (1 - xf 
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where 



and 



A u =2(C A i + C B i), 



A21 = CaiCbi, 



(16) 



(17) 



with Ca(b)i = cos (®a(b) +27rm;/3) and (m x ,m y ,m z ) = 
(1,-1,0). Of (I = x,y,z) is the relative angle between 
nearest neighboring spins in direction I and &a(b) * s 
the angle of the orbital pseudo-spin in sublattice A(B). 
(Tit) st is proportional to the band width W of the spin- 
and orbital-less fcrmion as 
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orders of M 2 and Mf as follows 



and is expanded with respect to M s and M t up to the 
lers o 
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x (l + a u M t + a M M? + a 3 /M t 3 ) , (19) 



N 



= - t ° x ^2 E (l + ^osO?M 2 

l=x,y,z 



where 



ot\i — - (Cai + Cbi), 



Q-21 



C 2 



Al 



cl 



BIJ 



4 

-CaiCbi, 



(20) 



(21) 



and 



4.33.4.3 2 
a u = iq$( C ai + c bi) ~ -^{CaiCbi + CaiC bi ) 

+ -L(Cai + Cbi). (22) 
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The detailed derivation of Eq. (|n]) is presented in Ap- 
pendix. It is worth to note that the terms which are 
proportional to M t or Mf appear in (Hj) st and (Tit) at- 
This is because the inversion symmetry in the system 
with respect to the orbital pseudo-spin is broken, i.e. the 



6 



free energy with T z — 1/2 is different from that with 
T z = —1/2. This is highly in contrast to the spin case 
where the inversion symmetry with respect to the spin 
operator is preserved due to the time reversal symmetry 
in the system without magnetic field. Being based on 
Eqs. (15) and (|l~9|), we investigate the phase transition in 
the F and A-type AF spin structures in more detail. 

Ferromagnetic structure: Eqs. ( |l5|) and (|l^) with the 
relation 9^ = for I = x,y and z are given by 



(Hj). 
N 



= - (1 - x) 2 -Ji(3 + M 2 )(l + B 2 M, 



- (1 - x)\ J 2 (l - M 2 )(3 - B 2 Mf), (23) 



and 



N 



= - tf)X- 



16 



x 1 + hMf + P 3 M, 



respectively, with 
B 2 = 



■ cos(9^ 



ei 



(24) 



(25) 



(26) 



in Sect. V. The first order transition in the cooperative 
Jahn- Teller system and its relation to the terms propor- 
tional to Q 3 , where Q indicates the normal mode of a 
MnC>6 octahedron, was discussed in Ref. [}?]. 

The parameter j3% becomes the largest at 9^ = 9 B = 
(2n + 1)tt/3 with n — (1,2,3) and determines the or- 
bital state uniquely. With taking into account j3 2 to- 
gether with /?3, (Tit) st favors the F-type orbital ordered 

(the so-called leaf- 



4 3.r 2 



and d 



3tf 2 



state with d x 2_ y 2, d y 2_ z 2 and d z 
type orbital) rather than d 3z 2_ 1 

(the so-called ciger-type one). Since the band width 
in the F-type orbital ordered state at T = docs 
not depend on types of the orbital, the thermal fluc- 
tuation stabilises the leaf-type orbital; the F-type or- 
bital ordered state with d x 2_ y 2 (c? 3z 2_ r 2) mixes with the 
AF-type one with {d 3x 2 _ r 2 / d 3y 2 _ r i) {{d y 2_ z 2 / d z 2_ x 2)) 
through the thermal fluctuation. The band width in the 
(d 3x 2_ r 2 /d 3y 2_ r 2) state is 5io which is larger than that in 
the [d y 2_ z 2 1 d z 2_ x 2~) one (3to)- The smaller band width in 
the latter state is attributed to the fact that the hopping 
integral in the xy plane is zero in this orbital state. 

A-type AF structure: In this spin structure, {Hj) s t 
and (Ht)st with Q s x = 9*=0 and 9* = n are given by 



N 



1 



(i- x y—(j 1 D 1 + j 2 D 2 ), 



(28) 



with 



and 



£>i = -18-2M. 



il^cos^ 
^r(cos9^ + 

^(cos(29^ 



9 f R 



cos 9, 



9 



cos(9^ + 29^ 



(27) 



The terms being proportional to M t vanish due to the 
cubic symmetry of the ferromagnetic spin structure. The 
coefficients B 2 and f3 2 become the largest at 9^ = 9 B — ir 
and at 9^ = 9 B , respectively. That is, Hj and TL t favor 
AF- and F-type orbital ordered states, respectively. 

Let us focus on the term being proportional to Mf in 
(7~Lt)st- In the case where this term is relevant, the orbital 
order-disorder transition is of the first order according to 
the Landau criterion in the phase transition. It corre- 
sponds to the transition at Too m the region of x > 0.25 
in Fig. 1 and is consistent with the first order transitions 
at the orbital ordering temperature observed in several 
manganites. On the other hand, the transition in the re- 
gion of x < 0.25 in Fig. 1 is of the second order. This is 
because 1) the term being proportional to M 3 does not 
appear in [Hj) s t, and 2) in this hole concentration re- 



ft' 



9* - 



gion, the G-type orbital ordered state with 9 
is realized. The inversion symmetry with respect to the 
pseudo-spin operator is preserved, i.e. /3 3 = in Eq. ( p7| ) 
in this orbital ordered state. The phase transition in 
this hole concentration region is discuss in more detail 



9cos(9^-9 B )Af 2 
(3 cos(9^ - 6 B ) - 2 cos 9^ cos 9 B ) 



M;m, (29) 



and 



D 2 = -18 + 6M 2 - 3 cos(9^ - 9 B )M 2 
- 4(cos9 A + cos9 B )M 2 M t 



and 

(Ht)s 
N 



with 



!3 cos(9^ - 9 B ) - 2 cos 9^ cos 9 B ) M'^M't, (30) 



16 



= -tnX\ —r 



x (l + 1m 2 + 7l M 2 M t + /3 2 M 2 + /3 3 M t 3 ) , (31) 



71 



= - — (cos 9*i + cos 9^ 
15 



(32) 



respectively. The most remarkable difference of the re- 
sults in the A-type AF structure from that in the ferro- 
magnetic one is the terms being proportional to M 2 M t . 
The origin of these terms is the anisotropic spin structure 
in the A-type AF state which breaks the cubic symmetry 
in the system. Because of these terms, the order param- 
eter of the A-type AF state acts as a magnetic field on 
the orbital pseudo-spin space and the relation Too > Tn 



7 



is derived. This relation is seen in the present phase di- 
agram (Fig. 1) and also in the experimental results in 
several manganites. 

Now we focus on these terms being proportional to 
M t Mg. At x — where Hj is dominant, the AF-type 
orbital ordered states with (O*/0* + 7r) for any 0* is real- 
ized above Tjv- Below Tjv, the term being proportional to 
M 2 M t becomes relevant and the orbital state is uniquely 
determined as (0*/ — O 4 ) with 



3(3 Ji - J 2 )M t + [J x + J 2 )M*M t v ; 

With decreasing temperature below Tjv, 0* continuously 
changes from tt/2 at Tjv to cos _1 {2 J 2 /(5 Ji — J2)} at 
T = 0. This orbital state favors the anitiferromagnetic 
interaction in the z direction in this spin structure. In 
the highly doped region where Ht is dominant, the term 
71 M s 2 M t in Eq. ( j3l|) favors the F-type orbital ordered 
state with Q l A = & B = n which corresponds to the state 
with d x 2_ y 2, as mentioned in the previous section. The 
terms being proportional to M^M t also appear in the 
C-type AF spin structure where the relation Too > Tjv 
is derived. The coefficients of the term in (TCj) s t and 
(Ht)st are given by J 2 {1 — x) 2 (cos0^ + cos0^)/2 and 
— t x32(cosQ* A + cos0*g)/(57r 2 ) which favor the G-type 
with (ds z 2 _ r 2 1 d x 2 _ y 2) and F-type with d 32 2_ r 2 orbital 
ordered states, respectively. 



region of x < 0.25 is of the second order as shown in 
Fig. 1. With taking into account the higher order cou- 
pling between the pseudo-spin and the Jahn- Teller type 
distortion in a MnOe octahedron and the anharmonic 
term of the potential energy for the lattice distortion, 
the phase transition changes from the second order to 
the first oneED However, since the phase transition ex- 
perimentally observed in LaMn03 at 780K is almost of 
the second order JIlMi 2 ] it is supporsed that the effects are 
small or are canceled out with each other in the com- 
pound. 2) The relation Tc > Too is satisfied in highly 
hole doped region (x > 0.1). This is because gain of 
the kinetic energy of electrons accompanied with the or- 
bital ordering is lower than that with the ferromagnetic 
one due to the hopping integral between different kinds 
of orbital unlike spin case. 3) In the region where Too 
and Tq are close with each other, the novel phase tran- 
sition is brought about due to the coupling between the 
spin and orbital degrees of freedom. The magnetization 
is enhanced below Too and Too increases by appling the 
magnetic field. These results well explain the unique ex- 
perimental results observed in Lao.8sSro.i2Mn03. 4) The 
sequential phase transition from the A-type AF phase to 
the ferromagnetic one with increasing T is caused by the 
thermal fluctuation of the orbital from the d x 2_ y 2 orbital 
ordered states. The ferromagnetic interaction in the z 
axis becomes finite due to the orbital fluctuation. 



V. SUMMARY AND DISCUSSION 

In this paper, we study roles of the orbital degree of 
freedom in phase transition in perovskite mangnanites. 
The effective Hamiltonian which includes the orbital de- 
gree of freedom as well as the spin and charge ones is 
utilized and the mean field theory at finite tempera- 
ture and carrier concentration is adopted. Through both 
the numerical and analytical calculations based on this 
theory, it is revealed that several characteristics of the 
phase transition observed in manganites originate from 
the unique properties of the orbital degree of freedom. 
The obtained results are summarized as follows: 1) The 
orbital order-disorder transition is of the first order in the 
wide region of x, and Tjv for the anisotropic spin struc- 
ture, such as the A- and C-type AF ones, is lower than 
Too- Both the results originate from the fact that the 
inversion symmetry in the system is broken with respect 
to the orbital pseudo-spin operator and the terms being 
proportional to M 2 M t and Mf exist in the free energy. 
These results are consistent with the phase transition oh-. 
served in Pr .5Sro.5Mn0 3 and Nd .45Sro.55Mn0 3 EME!l 
where the A-type AF state with d x 2_ y 2 orbital ordered 
one appears. The calculated results may be also appli- 
cable to the orbital ordering associated with the CE- 
type AF one which is experimentally confirmed to be 
of the first order transition in Pri_ x Ca x Mn03.E3 In the 
present calculation, the phase transition at Too in the 
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APPENDIX: EXPANSION OF (H t )st 

In this appendix, we present derivation of Eq. (JTsj) , 
i.e. the expansion of (H.t) s t with respect to the spin and 
orbital order parameters. We start from Eq. ( jlq ) where 
the band width W is calculated from Eq. mj) as follows 

w = 2j2r s s n, (Ai) 

s 

where 

ii = j2(K^\) s > ( A2 ) 



8 



and 



n = u 



if ij 3l'\lt 



(A3) 



Here, 5 indicates a vector connecting site i and it nearest 
neighboring site j. The spin part i| has the same form 
with Eq. (27) in Ref. 153 and is given by 



Ii = § cos Q s s Ml 



(A4) 



where 9| is the relative angle between spins at site i and 
site j and the relation M s = A s /3 + 0(A s3 ) is used. As for 
the orbital part, we present the derivation of I\ with S = 
±az termed I\ where a and z indicate a cell parameter 
of the cubic perovskite lattice and the unit vector in the 
z direction, respectively. 1\ with 5 = ±ax(y) is given by 
I* where 6* is replaced by 9* + 2tt/3 {Oj - 2tt/3). I* is 
calculated as 

9\ 9] 
cos — cos -f 
2 2 



v JO 



dSOi 



ej + so* 



d59 .g^*( COBI ^«+ cos <^j) 



(A5) 



where 



9\ —9*. The right hand side in Eq. (Ka) is 



expanded with respect to A* up to the order of 0(X CS ) as 
follows 



It 



to 



(2tt) 2 



1 



A t2 \2 



x ( Go + OiA 4 + ^A t2 ' Cj3 ^ t3 



2! 

where Q n (n — ~ 3 ) are given by 

Cio = 4, 



C»i = 3 cos 9*, 



3! 
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(8 - cos 2 9*), 



and 



-(cos 3 9* + 8 cos 9*). 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



The relation v l = 2tt(1 + ±A* 2 ) + 0(A* 4 ) is used. In the 
G-type orbital ordered state considered in Sect. Ill, 9' 
in the above formulas is replaced by 6^, B j, when site i 

belongs to the orbital sublattice A(B). B y u tilizing the 
relation M t = A*/2- A t3 /16 + 0(A* 5 ), Eq. @ is derived. 
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